www.gusucode.com > 网格法与lingo,内含例子及多个matlab代码 > 网格法与lingo,内含例子及多个matlab代码/third day/f3.m

    function y=f3(x1,x2,x3,x4,x5,x6,x7)
y=148257/1000*x1/x5/(x3/(x2-x1))^(3/20)*((1-131/50*(1-9/25/(x4/x2)^(14/25))^(3/2)*(x4/x2)^(29/25))/x6/x7)^(1/2)/(x2-x1);